2*2 fractionalized Chern/topological insulators in an exactly soluble correlated model 
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In this paper we propose an exactly soluble model in two-dimensional honeycomb lattice, from 
which two phases are found. One is the usual Chern/topological insulating state and the other is an 
interesting Z2 fractionalized Chern/topological insulator. While their bulk properties are similar, 
the edge-states of physical electrons are quite different. The single electron excitation of the former 
shows a free particle-like behavior while the latter one is gapped, which provides a definite signature 
to identify the fractionalized states. The transition between these two phases is found to fall into the 
3D Ising universal class. Significantly, near the quantum transition point the physical electron in 
the edge-states shows strong Luttinger liquid behavior. An extension to the interesting case of the 
square lattice is also made. In addition, we also discuss some relationship between our exactly soluble 
model and various Hubbard-like models existing in the literature. The essential difference between 
the proposed Z2 fractionalized Chern insulator and the hotly pursued fractional Chern insulator is 
also pointed out. The present work may be helpful for further study on the fractionalized insulating 
phase and related novel correlated quantum phases. 



I. INTRODUCTION 

The topological insulators/superconductors discovered 
recently, which cannot be described by the classic Lan- 
dau symmetry-breaking theoryji"— have motivated many 
studies in these novel states of matter liriSi The topologi- 
cal insulators are gapped in the bulk but have gapless he- 
lical edge-state protected by the time-reversal symmetry. 
Recently, the idea of topological insulators has been ex- 
panded into the symmetry-protected topological (SPT) 
states, which are bulk-gapped quantum phases with sym- 
metries, and have gapless or degenerate boundary states 
as long as the symmetries are not broken<^"— 

Another kind of state of matter beyond the Landau 
symmetry-breaking paradigm is the well-known integer 
quantum Hall (IQH) states and the fractional quantum 
Hall (FQH) statesii The former exhibits integer quan- 
tized Hall conductance in the external magnetic field 
while the latter shows the topological order and has ex- 
otic fractional charge excitations and fractional statis- 
tics. The Chern insulator is a band insulator exhibit- 
ing a nonzero quantized Hall conductance but preserv- 
ing the lattice translational symmetry. It is a natural 
extension of IQH in the lattice systems without any ex- 
ternal magnetic field. '^-'^'•^^ When including strong inter- 
action with appropriate electron-filling, fractional Chern 
insulator may exist in certain models and their nature 
has been extensively exploredi^"— 

In most cases, interaction effect in above mentioned 
topological states is not easy to treat, thus many ap- 
proximate analytical and/or numerical techniques and 
phenomenological approaches have to be utilized to ac- 
quire intuitive physics. Here we show a kind of exactly 
soluble models which may realize certain fractionalized 
Chern/topological insulators-like states in suitable pa- 
rameter space. Our proposed model is motivated by the 
study of the orthogonal metalS)^""— which have the same 



thermal and transport behaviors as the usual Landau 
Fermi liquid but with gapped single particle spectrum. 
Authors in Refs.|45l-l48| mostly focused on the issue of the 
exotic metallic states with the aim at possible non-Fermi 
liquid behaviors and the elusive critical Fermi surfaces 4^ 
Here we are interested in possible Z2 fractionalized states, 
which are hard to be identified by their bulk properties 
but their edge-states of physical electrons are gapped in 
contrast to their non-fractionalized counterparts, which 
can be used to identify these Z2 fractionalized states. 
This is an interesting feature shown in the present pa- 
per, which has not been reported in the previous works. 

The model we proposed is defined on various two di- 
mensional lattices at half-filling with the superficial Z2 
gauge structure. We provide an explicit formalism for 
the case of the honeycomb lattice and extend it to the 
interesting square lattice case. Using the dual trans- 
formation for the Z2 lattice gauge fields, the original 
model can be written as two decoupling Hamiltonian. 
One part describes a quantum transverse Ising model 
while the other is just free auxiliary fermion similar to 
the celebrated Haldane or Kane-Mele modelj^i^ The free 
auxiliary fermion always contributes quantized charge 
or spin Hall conductance while the dual quantum Ising 
model supports a second-order quantum phase transi- 
tion. The behavior of physical electrons is determined 
by combining these two parts and we find that the bulk 
properties are difficult to use for distinguishing Z2 frac- 
tionalized Chern/topological insulators from the usual 
Chern/topological insulating states. 

Significantly, the edge-states of physical electrons have 
rather different behavior for the fractionalized and the 
usual non-fractionalized states, which provides a defi- 
nite signature to identify the fractionalized states from 
the non-fractionalized ones. In the fractionalized states, 
since physical electrons lose their coherence, the cor- 
responding single electron excitation in edge-states are 



gapped while gapless edge-states survive in the usual 
Chern/topological insulating states because the auxiliary 
free fermions have non-zero weight of the physical elec- 
trons. (We should emphasize that although the single 
electron excitation is gapped in the fractionalized states, 
the static and dynamical many-particle correlations, e.g. 
density-density or spin-spin correlation, are still gapless 
in those states as what have been shown in Ref.|17|.^ 
More interesting, when the dual Ising model approaches 
its quantum critical point, the Green's function of phys- 
ical electron in the edge-states will show the strong Lut- 
tinger liquid behavior in contrast to the case of usual free 
particles. 



In addition, we have further made a comparable study 
to the fractionalized quantum spin Hall (QSH*) state in 
RefllTI and have extended our discussion to the case of 
the square lattice with a modified model. Particularly, 
for the square lattice, we expect that a chiral topological 
superconducting phase and its Z2 fractionalized version 
may appear when the attractive local interaction is intro- 
duced. Since our exactly soluble models may not appear 
in many real models for condensed matter physics, in 
term of the Z2 slave-spin mean-field approximation ; ^ ^ 1 ^^ 
we have suggested that similar fractionalized states may 
appear in various Hubbard-like models without intrin- 
sic degree of freedom for gauge fields. The relation of 
Z2 fractionalized Chern insulators found in the present 
work to the fractional Chern insulator Refs. [33 - l3a | is ex- 
plored, which shows that such two kinds of fractionalized 
states are rather distinct in their nature and construct- 
ing exactly soluble models for the desirable fractionalized 
Chern insulators is still widely open.— We also provide a 
brief discussion on the relation to orthogonal metals^"— 
and fractional topological insulators J^ 



The remainder of this paper is organized as follows. In 
Sec. nil we first introduce the basic exactly soluble model 
without spin degree of freedom on the honeycomb lattice 
and discuss two useful limit cases for this model. Then 
Sec. mil provides the solution of the basic model in terms 
of a dual transformation for the Z-i gauge-field. States 
with and without Zi fractionalization are found and one 
of them is identified as the Z-i fractionalized Chern insu- 
lator state while the other is the usual Chern insulator. 
The bulk features are distinct in the topological level, 
which emphasizes the confinement and deconfinement of 
two states. Properties of the corresponding edge-states 
are studied in detail. In Sec. IIVI the spin degree of 
freedom is included and the resulting states are also an- 
alyzed. In Sec. |Vl we present some possible extensions 
including the interesting soluble model on the square lat- 
tice and discussions on related issues. Finally, Sec. IVII is 
devoted to a brief conclusion. 



II. AN EXACTLY SOLUBLE MODEL 

We propose following model defined on the honeycomb 
lattice at half-filling, 

a = Hi + He, 

H. = -tY,4^t,c,~t!Y. 

(rj) ((ij)) 






(1) 



where Hj describes a modified Z2 lattice gauge theory 
model and He denotes the coupling of conduction elec- 
trons Ci to the former Z2 field af (Ising field). In Hj, h- 
term can be considered as the kinetic energy while J-term 
acts like a potential {j — i + 6a denotes three nearest- 
neighbor sites). [Readers who are not familiar with the 
lattice gauge theory may find Ref.jSJ] readable and use- 
ful.] He will be the standard Haldane model defined 
on the honeycomb lattice if Z2 gauge field is ignored. 
The first term of He is the usual hopping term between 
nearest-neighbor sites with tuning by the Z2 field a^,. 
The phase cpij = i^Tr in the next- nearest- neighbor hop- 
ping term in He is introduced to give rise to a quantized 
Hall conductance without external magnetic fields and 
the positive phase is gained with anticlockwise hopping. 
Here two Z2 gauge fields are introduced to enforce the 
Z2 gauge invariance and the label I in afjirf should be 
taken as the intermediate site between i and j sites. 

One may find our model is similar to the exactly sol- 
uble models provided in the study of the orthogonal 
metals4^"— Actually, this model is motivated by those 
works but we here focus on different respects. We are 
interested in the possible unusual gapped states while 
those authors mostly studied the metallic states in or- 
der to find possible non-Fermi liquid behaviors and the 
exotic critical Fermi surfaces. Since this model is insu- 
lating, its ground-state is hard to be identified by their 
bulk properties but one will find that the edge-states of 
physical electrons play a crucial role. This is not stud- 
ied in the previous works for the orthogonal metals and 
is an interesting feature in the present paper. In addi- 
tion, the authors in Ref.[13| also studied the unusual Z2 
gapped states (QSH*) with the Z2 slave-spin mean-field 
approach, we will discuss this point in Sec. IIVI and IVl 

Below we first analyze two simplified cases, which re- 
cover two models extensively discussed in the literature. 



A. The modified Z2 gauge theory model on the 
honeycomb lattice 

We first consider the first part Hi decoupled with the 
conduction electrons, the resultant Hamiltonian reads 



Hi^-hY^a:^ 



J 



E n 



(2) 



This model has the Z^ gauge symmetry (structure) and 
can be seen as fohows. First, one can define the Zi 
gauge transformation operator as Gi = Iljeheajagon ^Ij 
and this operator flips all Ising gauge field a^ in the i- 
th hexagon. Since J-term covers two Ising gauge fields 
in the same hexagon, this term is obviously unchanged 
when the Zi gauge transformation Gi is utilized. There- 
fore, the Hamiltonian ^ is invariant under the operation 
of Gi, thus it has the expected Zi gauge symmetry (in- 
variance). One can also check that \Gi^ili\ = 0, which 
means that the Hamiltonian has the local (gauge) sym- 
metry enforced by Gi. Since Gf = 1, one may choose 
Gi = ±1. Moreover, the physical states of the Hamilto- 
nian should be invariant under the operation of Gi . (For 
a brief discussion on the standard Z2 gauge theory model 
one can refer to Appendix A.) 

There are two limit situations which are easy to an- 
alyze. First, ii J ^ h, the ground-state should have 
n7=i-i-<5„ '^ij — 1 ^i^d this is just the deconfined state 
for the Z2 gauge theory. However, when J ^ h, the 
ground-state can be obtained by setting all (xf ■ = 1 and 
the resultant state is the confined state for the Z2 gauge 
theory. More importantly, any physical excitation must 
carry zero gauge charge in the confined state in contrast 
to the case of deconfined state where fractionalized exci- 
tation with nonzero gauge charge are permitted in prin- 
ciple. Such feature for the excitation will be manifested 
when matter fields (electrons) are included in the gauge 
theory model. 

To sharpen our understanding of the above model, it 
is useful to carry out a dual transformation which trans- 
forms the Z2 lattice gauge theory model into a simple 
(quantum) Ising model living on the honeycomb lattice. 
The dual transformation is defined as 



n 



Ti ,cr 



T- T ■ 

I] 'l '3 



Then, the dual Hamiltonian reads as 



H, 



(ij) 



jT.-' 



(3) 



which is just the well-known quantum transverse Ising 
model on the honeycomb lattice:- We should empha- 
size that the Z2 gauge transformation operator Gi = 
n^-i^x CTf, = 1 due to above dual transformation and 
no constraints are needed for the dual Hamiltonian in 
contrast to the original one. The ground-state phase di- 
agram of the quantum transverse Ising model is clear: it 
has two phases characterized by the usual Landau local 
order parameter (ff ) = (the paramagnetic phase) and 
(ff ) 7^ (the ferromagnetic phase) J' Besides, a second- 
order quantum phase transition exists between those two 
phases whose critical behaviors belong to the 3D Ising 
universal classj^ 

Some readers may wonder why we do not use the stan- 
dard Z2 gauge theory model on the honeycomb lattice, 
which is more friendly and is more easy to analyze its 



topological order (For a detail see Appendix A). The rea- 
son is that in the standard model, the Z2 gauge transfor- 
mation operator Gi cannot automatically be unit after 
the dual transformation and this leads to unnecessary 
extra complications. Therefore, we use the modified Z2 
gauge theory model which reflects the same physics with- 
out caring about the unsatisfactory constraints. How- 
ever, the modified Z2 gauge theory model is not friendly 
if one considers its topological order although the modi- 
fied Z2 gauge theory model itself indeed shows the same 
topological order as the standard one. 

Additionally, if one is only interested in the low-energy 
physics, the following ip'^ theory may be useful in this 
respect. 



Z 



D(t)e 



- J dTd^x[{d-r4>f +c^ (V cjif +r4,'^ +u4>^] 



(4) 



where r, u are effective parameters depending on micro- 
scopic details. (For the derivation and brief discussion 
on this effective action one can refer to Appendix B.) 

According to the above dual transformation, it is read- 
ily to see that the paramagnetic state of the quantum 
Ising model corresponds to the deconfined state of the 
original Z2 gauge theory model while the ferromagnetic 
phase dualizes to the confined state. Moreover, the quan- 
tum critical point of the quantum Ising model trans- 
lates to the deconfinement-confinement transition point 
though no Landau local order parameter can be reliably 
defined in the case of gauge theory model. (The Wil- 
son loop can be defined in the gauge theory model but it 
is not a local operator and we will not study this point 
further in our present paper.) 



B. The pure Haldane model 

In this subsection, we will briefly review some impor- 
tant properties of the pure Haldane model (without any 
gauge fields coupled to the electrons), which is useful for 
our next discussion. 

Our starting point is the usual Haldane model on the 
honeycomb lattice at half-filling'^^ 



H 



H 



-*E4c,-t'E 



.iVijJ 



C Ci 



(5) 



(y-> 



{{ij)) 



It is useful to rewrite this single-particle Hamiltonian in 
the momentum space as 

Hh = -tJ2[f{k)cl^CkB + nk)cUckA] 

k 

+ 2t'j{k){cl^CkA - clBCks), 

where we have defined f{k) — e^^'^^ + 2e'''^/^ cos(^fcy), 

^(A:) = sin(\/3fcy) — 2cos(|fca;)sin(-^A:j^) and A, B repre- 
senting two nonequivalent sublattices of the honeycomb 
lattice, respectively. Then, by diagonalizing the above 



Hamiltonian, one obtains the quasiparticle energy band 
as 

Ek± = ±Vt2|/(fc)|2+4i'27(fc)2, 

which preserves the particle-hole symmetry. It is well- 
known that for 3-\/3i' < t, the excitation gap mainly 
opens near six Dirac points (Only two of them are 
nonequivalent in fact):^^ Then, expanding both /(fc) 
and 7(fc) near two nonequivalent Dirac points ±K — 
±(0, Y%)^ respectively, the gap can be found as Agap = 

6v3t' and the quasiparticle energy reads Eq„± ~ 

±y^(|ig)2 -t- (3V3t')' with q = {q,,qy) = (fc„ fc, T ^)- 
The most interesting property of the Haldane model is 
that it gives rise to a quantized Hall conductance with- 
out applying the external magnetic field in contrast to the 
usual quantum Hall effect. The corresponding quantized 

2 

Hall conductance can be readily calculated by cth = Cij^ 
where Ci = ^ / dk^dkyd . (^ x ^),S For the Haldane 
model, the single particle Hamiltonian can be rewritten 
a.sh{k) ==d(fc)-o-withd = {~tRcf{k),tlmf{k),2t'j{k)), 
d = d/|d| and a being the usual Pauli matrices. Then, 
by inserting the expression of d into the formula of Ci, 

2 

one obtains Ci = 1 and an = 7;- for the Haldane model. 
In literature, the ground-state of the Haldane model is 
called the Chern insulator since it exhibits a quantized 
Hall conductance without external magnetic fields and 
preserves the lattice translational symmetry of the hon- 
eycomb lattice."--'"^ Hereafter, we will use the Chern in- 
sulator when involving the ground-state for the Haldane 
model. 



The low-energy effective theory for the pure 
Haldane model 



For the discussion of the low-energy physics, an effec- 
tive 2 + ID massive Dirac action can be obtained by 
expanding original Haldane model around two nonequiv- 
alent Dirac points ±K, 



coupled to the external electromagnetic field reads 



S 



H 



/ (fxdrCo = / 



d'^xdT^[ipaht^df_, 



■w)V'a], 



where 7^ = (Tz,T^,Ty) and d^ = {dr,dx,dy) with 
TzTTx,Ty the Pauli matrices. Here the same indices 
mean summation. We have introduced the effective mass 
771 — ~3\/3t' of Dirac fermions and set the effective 
Fermi velocity vp = |t to unit. The Dirac fields are 

defined as Ipia = [ciAa^CiBaV, '02<t = {c2Aa,-C2Bc7V 

and i/jaa = i'aalo with a = 1,2 denoting the states near 
the two nonequivalent Dirac points ±K = ±(0, -^y=) and 
T implying the transposition manipulation. 

If the external electromagnetic field A^, — {i(p, Aj.,Ay) 
is introduced by the conventional minimal coupling 
(9^ -^ dfj, — ieA^), the resulting effective Dirac action 



5 = / d^xdr^l'tpai'^^id^ - ieA^ 



)+m)i;a]- (6) 



By integrating out the Dirac fields, we get an effec- 
tive Chern-Simons action, which represents the nontrivial 
electromagnetic response of the massive Dirac fermions 
to the external electromagnetic field y4^^(For details, 
see Appendix C.) 

Scs^ [d^xdT[Ne^^^e''''^Af,d,A^], (7) 

J 8n\m\ 

where N = 2 (from the two nonequivalent Dirac points) , 
g^i^A jg ^Yie usual all- antisymmetric tensor and the regu- 
lar Maxwell term (~ F^,^) has been dropped out since 
the low energy physics is dominated by the Chern- 
Simons term alone. Then, the physically observable 
quantized Hall conductance can be obtained from J^ — 



as. 



£s| ^ _ njp2 



TH^iA^o - ^-^ 4^(^y^o - d^Ay) = ^Ey and we 

2 
reproduce the result for the Hall conductance cr/f = ^ 

where h — 2T:h is reintroduced with 771/ 1 777 1 = — 1 and 
N ^2A 

Alternatively, one may inspect the low-energy physics 
from the perspective of the corresponding edge-state. 
Firstly, the following effective Chern-Simon action re- 
produces the quantized Hall conductance of the Haldane 
model 

S= f d^xdt[^e^''^a^d,ax + ^ef^-'^A^d^ax], (8) 

J 47r ZTT 

where a^ is the auxiliary dynamic gauge field and it 
gives rise to the physical charge/particle current j^ = 
i^t^^^d^ax. According to the formalism in quantum 
Hall effectfi the so-called K-matrix in the present ac- 
tion \s K = 1 with the charge- vector g = 1, thus 

the physical quantized Hall conductance is calculated as 

1 2 2 

an = qK" q^ = 77 ■ Moreover, the ground-state of 

the present model is not degenerated in the torus due 
to \K\ — 1 and the elementary quasiparticle is the usual 
fermion (electron) since the the exchange statistical an- 
gle is TT, which means that the exchange of two iden- 
tical quasiparticle leads to a tt phase acquired in their 
wavefunction. The corresponding edge-state can be eas- 
ily derived by standard bulk-edge correspondence for the 
effective abelian Chern-Simons theory," 



S. 



edge~ 



dxdt-—[dt4'dx(j) - v{dx(i)f] 
47r 



---{dt(j)Ax -dx(i 
III 



(9) 



with V denoting the non-universal velocity of edge states 
and (/> being the bosonic representation for the edge-state 
modes. This edge-state can be refermionized by introduc- 
ing the fermion operator ij} oc e**^ and the resulting action 
reads 

Sedge = / dxdt[ip\idt + eAt - ivdx - veAx)^!)]. (10) 



The above edge-state mode contributes e^ /h to the quan- 
tized Hall conductance and its gaplessness or the exactly 
quantized value for the Hall conductance is protected by 
the chiral feature of the edge-state. 



III. THE CHERN INSULATOR AND THE Z2 
FRACTIONALIZED CHERN INSULATOR 

Having studied the pure Haldane Eq.(l5]) and the Z2 
gauge theory model Eq.Q, it is ready to discussion the 
full model Eq.([T|) with both gauge field (Ising field) and 
matter field (electrons). For convenience of discussion, 
here we rewrite Eq. ([1]) as follows 

(ij) i j=i+Sa 

{«> ((y» 

Comparing this model to the Z2 gauge theory model 
Eq.® and its related Z2 gauge transformation operator, 
one may use the same Q = Ilje/ieajagon '^Ij^ which leads 
to the happy result [Gi,H] = 0. Thus, the full model 
Eq.(IT]) have the desirable Z^ gauge structure. One may 
note that the physical electrons do not carry any gauge 
charge of the Ising field af , which means even though the 
Ising field af. appears in its confined state, the electrons 
are still well-defined physical excitations in this case. 

Here, let us first inspect the model ([T]) without resort- 
ing to its dual formalism. Since the physical electrons are 
gapped due to the next-nearest-neighbor hopping, one 
may integrate them out and the resulting Hamiltonian 
will be like the effective pure gauge field part^l 



{ij) 



^E n ^i 



V 



n 



^-^ J. J. ' V 

j^hexagon 



If the Z2 gauge-field is in its deconfined state and V is 
small, one may set all ni=i;+a„ ^fj — 1 ^^id neglect other 
two terms. Then the low-energy excitation in this case is 
the so-called Z2 charge ^ii The Z2 charge near site i can 
be seen as a manifold with TTo-.j-a o'?:, = —1 with others 
unchanged. When V becomes large, the system remains 
in the deconfined phase but there exists another low-lying 
excitation the Z2 vortex (vison), which is created by set- 
ting certain Yljahexagon^h = -1 while keeping others 
intact*^ It is easy to see that /i-term adds or destroys 
the Z2 charge while V-ieiuY do not change the number 
of Z2 charge. In addition, the deconfined Z2 gauge the- 
ory has the so-called topological order, (A characteristic 
signature of the topological order is the ground-state de- 
generacy depending on the topology of the system) e.g., 
there exists four (two) degenerate ground-states if the 
system is put on a torus (cylinder) 1^ 



We also note that an effective mutual f/(l) x f/(l) 
Chern-Simons theory can describe the low-energy behav- 
ior of the Z2 gauge theory modeli^ 



S = d xdr 



Since the /C-matrix is K 



and iDetKl 



-e'^'^a^d^bx + a^j^, + 6^ J^ 



where j^ and J^ represent the Z2 charge and the Z2 vor- 
tex, respectively, a^, 6^ are the auxiliary gauge fields 
whose Chern-Simons term represents the semionic mu- 
tual statistics between the Z2 charge and the Z2 vortex. 

2 

2 

the mutual U{1) x C/(l) Chern-Simons theory correctly 
reproduces the four degenerate ground-state on a torus. 
When the gauge-field is confined, the /i-term will dom- 
inate and the number of Z2 charge is not well-defined. 
In other words, the Z2 charge is condensed in this situ- 
ation and we may consider the confined state is just the 
condensed state for the Z2 charge. The ground-state in 
the large h limit is obtained by setting all fjf- = 1. In- 
stead, the elementary excitation is the Z2 link created by 
(jf^ = —1. In addition, the previous mutual C/(l) x U{\) 
Chern-Simons theory is not useful in this state due to the 
confined nature. 

Then, we will immerse into the discuss with the help of 
the dual transformation. Following the same treatment 
in the last section, the dual transformation is defined by 
nj=i-i-5„ '^fj = ^i^^ij = ^i^f ^-iid the resulting Hamilto- 



reads 
H = Hi + He, 

^ = -iE4^^c,-i'5]e'^"^clffffc,. (11) 

The above Hamiltonian can be further transformed via 
fi = ff Ci, {—lyi'^^T^ — > ff and one arrives at a simple 
formalism 

H = Ht + Hf, 

Hf = -tJ2flf,-fY.^''''f'fr (12) 

(■y) {{ij}) 

Obviously, Hf is identical to the usual Haldane model 

([5]) except that the / fermion is not the original electron 

c but acts like a slave-particle. Hr is just Eq.([3]) obtained 

in Sec. [Hi Since the Haldane model has a quantized Hall 

2 
conductance, Hf also has such Hall conductance (Jh — 7^ 

and the edge-state is described by Eq. ^ or Eq. (ITUl) . The 

topological properties of Hf will be reliably captured by 

Eq.© or Eq.®, as well. 



A. The nature of two insulating states 

Since the / fermion is always insulating and has the 
quantized Hall conductance, the ground-state of the 
whole system depends on the Ising field part Eq.(j31). 
From the discussion in previous section, we know the 
Ising field part can appear in the ferromagnetic, param- 
agnetic and quantum critical state. In the case of fer- 
romagnetic state (f^) 7^ 0, a ~ {T^)fi and the physical 
electron behaves as the / fermion, which can be identi- 
fied as the usual Chern insulator as what has been done 
in Sec. [Ill However, if the Ising field is paramagnetic 
((f^) = 0), the physical electrons will not be a well- 
defined quasiparticle but the / fermion and Ising field 
will be useful quasiparticle excitations. In this case, the 
Z2 gauge structure is meaningful since the elementary ex- 
citations are the / fermion and the Ising field, which both 
carry the Z2 gauge charge. Therefore, we may identify 
this state as a Z2 fractionalized state. When consider- 
ing that such Z2 fractionalized state is still insulating 
and has the quantized Hall conductance, we may call it 
the Z2 fractionalized Chern insulator. We note that the 
paramagnetic state corresponds to the deconfined state 
of the original gauge field and in this case fractional- 
ized excitations, e.g. / fermion and f^, are meaningful 
and can be choose as the real quasiparticle. Instead, 
the ferromagnetic state is dual to the confined state of 
the Z2 gauge theory and only physical electron c itself 
is the real quasiparticle excitations. Therefore, the dual 
Hamiltonian can reproduce correct results as expected 
from the original model. When the Ising field part is 
quantum critical, the system will approach its quantum 
transition point and the critical behaviors will be deter- 
mined by the Ising field part of Eq. Q, which means 
that the transition from the usual Chern insulator to the 
Z2 fractionalized Chern insulator falls into the classic 3D 
Ising universal class. (The critical exponents for the 3D 
Ising universal class are a = 0.11, /3 = 0.32, 7 = 1.24, 6 = 
4.9,1^ = 0.63,77 = 0.04,2 = 1.) 

It is interesting to see that the two insulating states 
(the Chern insulator and the Z2 fractionalized Chern 
insulator) discussed above are rather different but they 
cannot be distinguished by Landau symmetry-breaking 
theory and the usual topological quantum number such 
as the (first) Chern number. Therefore, one has to find 
other methods to identified these two phases when more 
general and complex models are encountered. A straight- 
forward idea is to check the Green's function of physical 
electrons in bulk but it does not work since they are both 
insulating and their Green's function of physical electrons 
are both gapped. However, the single particle Green's 
function of physical electrons in the edge may solve this 
issue since the Chern insulator will have a gapless edge- 
state while the Z2 fractionalized Chern insulator only has 
a gapped one. We will visit this important issue in the 
following subsection. 

Before proceeding, we should emphasize that the 
charge of original electrons c is only carried by the / 



fermion since f^ is a real operator and cannot carry the 

electromagnetic C/(l) charge. Therefore, the quantized 

2 
Hall conductance an — ^ obtained from Hf is the true 

one. 



B. The edge-state of the Chern insulator and the 
Z2 fractionalized Chern insulator 

The single particle Green's function for the physical 
electron is defined by G{x,t) = —i{Tc{x,t)c^ {0,0)) — 
-i(T/(x,i)/t(0,0)f^(x,t)f^(0,0)). Since the / fermion 
and Ising field f^ are decoupled as can be seen in Eq. ([T2|) , 
onehasG(a;,i) = -i(T/(a;, t)/^(0, 0))(Tf^(a;,i)f^(0, 0)). 
Recalling that the edge-state for / fermion is described 
by Eq. pUj) while the low-energy effective description of 
the Ising fields is the (^'* theory, the Green's function 
for the physical electron in the edge will behave as 
G{x,t) ex j^^ X G^. G^ is the Green's function of 
1+lD (y9^ theory which may describe the Ising field ex- 
citation near the edge. When the system is in the Z2 
fractionalized Chern insulating state (the Ising field is in 
its paramagnetic phase), G^ should show a gap, which 
leads to a gap for physical electrons and G{x, t) decays 
out for long time and distance. In contrast, when one 
considers the usual Chern insulating state, the Ising field 
is ferromagnetically ordered in this case, thus we may 
have G^ ex constant, which means that the physical 
electrons will have the same behaviors as the / fermion 

(G(x,t)cx^). 

Interestingly, if the Ising field is critical (this case is re- 
alized when the bulk is at its quantum transition point.), 
the Green's function of physical electrons in the edge will 
behave as 



G{x, t) ex 



x-vt [a;2 - (rf)2]'(/2 



(13) 



with G^ c>c [x"^ — [ct)"^]^^ . [rj = 1/4 denotes the anoma- 
lous dimension of the 2D classical Ising model and c is 
the edge- velocity for the Ising field.) In the long distance 
limit, we have G{x) ex x"^'"*, which decays more rapidly 
than the free electron G{x) ex a;~^. The corresponding 
local density of states is 7V(w) ex |w|^/^. In some sense, 
such explicit Luttinger liquid behaviors, which could be 
measured by local differential conductance in scanning 
tunneling microscopy (STM), may be used to locate the 
exact position of quantum critical point of the system, 
although the bulk critical fluctuations will dominate in 
the case. 

We should emphasize that although the single electron 
excitation is gapped in the edge-state of Z2 fractionalized 
Chern insulator, the static and dynamical density-density 
correlations are still gapless in those edge-states as what 
have been shown in Ref . [171 ^ This point can be seen 
as follows. The physical static density-density (charge) 
correlation function in the mean-field approximation 
reads (ct(x)c(a:)ct(0)c(0)) « (/t(x)/(x)/t(0)/(0)). 



Since / fermion is gapless in the edge, the static 
density-density correlation is expected to be gapless. 
For models including the spin degree of freedom, 
(Those models will be presented in next section.) the 
static spin-spin correlation will also show the gap- 
less behavior in the edge-state of Z2 fractionalized 
states. For the case of the dynamical two-particle 
correlation function, e.g. the dynamical density- 
density correlation (c'''(x, t)c(a;,i)c^'(0, 0)c(0, 0)) sa 
(/t(2;,t)/(2;,t)/t(0,0)/(0,0)), should also be gapless. 

After all, the edge state of physical electrons clearly 
show distinct behaviors in the Chern insulating and the 
Z2 fractionalized Chern insulating states. Therefore, one 
can check the behaviors of the single electron excitation 
in edge-state to identify those two insulating states in 
generical models. 



IV. EXACTLY SOLUBLE MODELS WITH SPIN 
DEGREE OF FREEDOM 

A. Extension to include spin degree of freedom 

Since most real models have the spin degree of freedom, 
it is crucial to include this element in our model. A 
simplest extension is to double the Haldane model part, 
which means electrons carrying different spin-flavor are 
not interacting at all. 

a = Hj + He , 

(ij) i j=i+Sa 



(ii>o- 



{{^j))'y 



Since electrons with different spins do not interact, prop- 
erties of this model is identical to Eq. ((TJ expect for 
doubled quantized Hall conductance and two chiral edge- 
states. For completeness, we also write down the dual 
Hamiltonian, which can be readily derived in terms of 
the same treatment in last section. 

H = Hr + Hf, 

(ij) i 

Hf = -tJ2 flf^'^ - ^' E ^""flf^^- (15) 



B. Z2 fractionalized topological insulator 

The next extension for including the spin degree of 
freedom is to replace the Haldane model part by the well- 



known Kane-Mele mo del ^ 
H = Hi + He , 

(ij) i j=i+Sa 



{ij)c 



{{^j))<y 



where the next-nearest neighbor hopping term is spin- 
dependent, which mimics the spin-orbit coupling effect in 
the real materials. The corresponding dual Hamiltonian 
reads 

H = Hr+Hf, 

(ij) i 

^f--^H flfi- - ^' E ^e'^'^/l/.- (17) 

{ij}a- ((«i)>o- 

For this model, we expect that the system can show 
two distinct phases, one is the usual topological insulator 
(quantum spin Hall insulator) ((ff) 7^ 0) and the other 
is the Z2 fractionalized topological insulator ((ff ) = 0). 
These two phases will have no quantized charge Hall con- 
ductance [an = cTff + (Tff = 0) but show quantized spin 

Hall conductance (cfsh — ^i'^H~'^H) — ^)- For the Z2 
fractionalized topological insulator, its low-energy prop- 
erties are described by the same pure gauge-field theory 
as what has been discussed in previous section 



H 



/--/^E-^-^E n ^^+^E n < 



iv) 



i j=i+5a 



i j£hexagon 



For their edge-states, the usual topological insu- 
lating states has the helical edge-state Hedge = 

J dx['ipl{—ivdx)4'^ — 4'l{—'i'V9x)ipi] and this gapless edge- 
state is protected by the time-reversal symmetry. For 
the case of the Z2 fractionalized topological insulator, 
the single particle excitation of physical electrons in the 
edge-state is gapped due to the massive Ising spin field 
f^. 

Besides, when the system is critical, the edge-state 
for the physical electrons has the behavior described by 
Ea. (ll3p with two spin-fiavors. One may write down the 
corresponding Luttinger liquid Hamiltonian for such crit- 
ical case^^ 



He 



1 

2^ 



dx[uK{da:0)^ 



K 



{d.c^f 



where u denotes the non-universal velocity and K is 
the Luttinger parameter. The fermion (electron) oper- 
ator is defined by tpcrix) oc e"'*^'^^^^"'^^^^-'^. One may 
find that the Luttinger parameter AT = 2 or AT = 1/2 
both can reproduce power-law behavior of the electron 
Green's function. However, generally, AT > 1 means the 
electrons are very likely to be pairing and considering 
no such instability exists in the original model, we con- 
clude one should choose K — 1/2 and the local den- 
sity of state for the physical electrons is N{uj) oc |w|^/^. 



The related momentum distribution function of physical 
electrons is n{k) ^constant+|fc|^/*s(7n(fc) and one clearly 
sees no jump appears near the Fermi momentum (fc = 0). 
The density-density correlation in this edge-state is easy 
to find and it decays as - in the large distance limit while 
the free electron case gives a -^ behavior. 

In Ref.[U, authors argued the existence of the frac- 
tional topological insulators, which can be realized by 
requiring that the two spin species in the usual topolog- 
ical insulator each form fractional quantum Hall states. 
The name "fractional" topological insulators results from 
the natural extension of the usual fractional quantum 
Hall effect and its corresponding topological order is ex- 
pected to be the similar one in the fractional quantum 
Hall states. For our Z2 fractionalized topological insula- 
tor, we have the Z2 topological order and do not involve 
the fractional quantum Hall-like objects. Thus, we think 
that the Z2 fractionalized topological insulator is differ- 
ent from the fractional topological insulators proposed in 
Ref.iH. 



C. Comparing the Z2 fractionalized topological 
insulator to the QSH* state 

In Ref . |17l| . an unusual Z2 gapped states (fractional- 
ized quantum spin Hall state (QSH*)) is discovered by 
using the Z2 slave-spin mean-field approximation on the 
honeycomb lattice. Such state is stable when local repul- 
sive energy and spin-orbit coupling are both strong and 
it has quantized spin Hall conductance as the Z2 frac- 
tionalized topological insulator studied in the last sub- 
section. In the low-energy limit, QSH* is described by a 
pure Z2 lattice gauge theory while the Z2 fractionalized 
topological insulator can also be described the pure Z2 
gauge theory. (One may simply discard the electrons c in 
Eq. ()15|) since they cannot contribute singular correction 
to the Z2 gauge fields.) Furthermore, the single electron 
excitations in edge-states for QSH* and the Z2 frac- 
tionalized topological insulator are gapped. Therefore, 
we may think the Z2 fractionalized topological insulator 
found here is similar to the QSH* state in Ref. [131 ■ I'^ 
some sense, our model may provide a possible realization 
for QSH* state in exactly soluble models. 



V. EXTENSIONS AND DISCUSSIONS 
A. Relation to the real microscopic models 

The models Eqs.p|).p4 |) . (fT6|) introduced in the present 
paper have the Z2 gauge field part, which is rather arti- 
ficial and dose not appear in many real models for con- 
densed matter physics. Therefore, it is important to find 
some real and simple models which can provide the sim- 
ilar states as what have been found in our main text. 

The first model is the topological Hubbard model on 



the honeycomb lattice at half-filling. 



55 



H=-tJ2i4.c,a + h.c.)^t' Y. e'^-^cL 



{ij)<^ 



{{^J))<^ 



u 



^E(-^-l)^ 



(18) 



where rii = X^o- "^L"^*'^' ^ ^^ ^^^ onsite Coulomb energy 
between electrons on the same site and t is the hopping 
energy between nearest-neighbor sites. Since we are in- 
terested in the case of half-filling, the chemical potential 
has been set to zero. 

Then in terms of Z2 slave-spin representation, ^^'^^ the 
physical electron Ca- is fractionalized into a new slave- 
fermion fa- and a slave-spin r^ as Cia — fiuTf with a 
constraint rf + 1 = 2{ni — 1)^ enforced in every site. 
Under this representation, the original Hamiltonian can 
be rewritten as 






J a 



-tE(-^ + i) 



(19) 



,/ - 



where m = nl ^ Y^a fLfi<y Obviously, a Z2 local 
gauge symmetry is left in this representation and the 
corresponding low-energy effective theory should respect 
this. The mentioned gauge structure can be seen if 
fll^ -)■ eifll^ and t[ -)■ tirf with e^ = ±1 while the 
whole Hamiltonian H is invariant under this Z2 gauge 
transformation. 

When one utilizes the mean-field approximation, the 
resulting mean-field Hamiltonian will be 

^/ = - E ^i^3flf3^ + h.c.) Y. ^l.e^'''^4c,., (20) 



(*jV 



((■y»'^ 



U 



Hi = -J: {J^^rrr^ + h.c.) J2 4<-f + T E«« 



{v)<7 



(fe)> 



f(^-T-) 



where we have defined i^ = t{TfT^),t[- 

Jn = tEMlM, J[, = t'e^^'^ZJflf.a). The de- 
coupled Hamiltonian Hi is an extended quantum Ising 
model in transverse field and Hf describes / fermions in 
the honeycomb lattice. Here, all the Lagrange multipliers 
have been set to zero, provided only non-magnetic solu- 
tions are involved and a half-filling case is consideredi^ 

Now, comparing Eqs. (pn)) . (PT]) to Eq. p5)) . it is obvious 
that the mean-field Hamiltonian is able to capture the 
basic features of the finely tuned exactly soluble models 
except for the extra trivial next-nearest-neighbor cou- 
pling term in the quantum Ising model part. However, 
one should realize that the mean-field treatment is a bi- 
ased method and the original model may not have to 
support the putative Z2 fractionalized states. 

If one wants to go beyond the mean-field approxima- 
tion, the intrinsic Z2 gauge structure should be respected 



and the above mean-field Hamiltonian could be modified 
into the following ones 

Hf = -J2{Ujflgtjfja + h.c.)- J2 *lje''^"4^';ft>j-> 

{ij)a {{ij))a- 

Hl^-J2 {J^3T^Q13t'^ + h.C.) - Y^ J[jTf QliQfjT^ 



Hz^ — 



i 

-hY^gf^-GY: n ^^ 



(22) 



(u> 



i j^hexagon 



where pf ■ denotes the dynamical Z2 gauge-field, which 
is introduced to enforce the desirable intrinsic Z2 gauge 
structure. For convenience, we have added the Z2 gauge- 
field part Hz2 and the gauge transformation operator 
is defined by G, = (-l)^r-HE„ //„A. ll,=^+s^ 9-y If 

we focus on the fractionalized states, the Hi could be 
dropped out and the gauge transformation operator is 

now read as Gj — (— l)^"-'^'"-'^*'' nj=i+5„ ^fj- Because 
Gf = 1 and [Gi,Hf + Hz^l = 0, we may set Gi = 1 for 
physical states. Then, it is straightforward to see that 
if a Z2 charge (n,=i+5„ Bfj = ~1) is created at site i, 
this site must be occupied by a slave-fermion f„ so as to 
fulfill the constraint Gi — 1. In contrast, if no Z2 charge 
appears, the corresponding site should not have a slave- 
fermion fa-. Therefore, we may say that the Z2 charge is 
bound to the slave-fermion and the mutual statistics of 
slave-fermion and the Z2 vortex iUj^^hexagon 0ij = "1) 
is the semionic statistics. (Recalling that the semionic 
mutual statistics between the Z2 charge and the Z2 vor- 
tex is discussed in Sec lIIII ) More detailed properties on 
similar model has been investigated in Ref.[17] and we 
refer interested reader to their original paper. 

In addition, when one wants to describe the Z2 frac- 
tionalized topological insulator in Eg. dTSl) . the simplest 
model will be given by 

H=-tY(4.Cja+h.C.)-t' Y <^^'^''4.C,a 

+ yE("^-l)'- (23) 



For the electron part, a simple and feasible choice is^ 
+H.C. + Y ^z{4^c^^ - 4^c^^) 

i 

+t'Y^4t^it-4i^3i) (25) 

where i — {ix,iy) and only nearest- neighbor- hopping is 
involved. rUz denotes a tunable parameter which may re- 
sult from the usual ferromagnetic order. We also require 
the implicit half-filling condition for this model. Follow- 
ing the treatment in the main text, the electron part will 
be dual to a free auxiliary fermions part, whose Hamil- 
tonian reads 

Hf^~YflMk)-gU,fk. 

k 

with gx, Qy, Qz being the usual Pauli matrix act- 
ing on the spin space. We have defined d(fc) — 
(2t sin fey, 2isinfc:c, — m^ -|-2i'cosfca; + 2i'cosfcy). As what 

has been done in the Haldane model. Hall conductance 

2 
of this model is calculated by an = C'lx with Gi = 

-J-Jdk^dkyd ■ (^ X §^). One finds Gi = 1^ for 

< |mz| < At' while Gi vanishes in other cases. Such 
Hall conductance can also be obtained by expanding the 
above Hamiltonian around four high-symmetry points of 
the energy band. ((0, 0), (tt, 0), (0, tt), (tt, tt)) 



S ■ 



(fixdrYVl'ailti.id^, 



ieA^) -I- ma)i}a] 



where the effective mass has been defined as to„-(o 0) 



At' — 7712 and lUa^ 

If 7712 — 0, TOq 



.(7r,0) 



4t'-|-TO2 when At' close to |7tt.2|- 
"^a=(o.7r) = ^712. We note the 
effective Dirac theory only has two-flavor, thus Gi can 
take only 0, ±1. Other properties are all similar to the 
case of the previous discussion and we will not present 
them here. 



C. The model for the topological superconductor 



B. The exactly soluble model on the square lattice 

In the square lattice, the gauge-field part is now read 
as 

Hi = -hY^'s ~ JT. n -S (24) 

(ij) ^ j—i±x,idzy 

where the Z2 gauge transformation operator is defined as 

^i ~ lljeplaquett '^ij' 



We have discussed the usual/fractionalized 
Chern/topological insulators in the previous sections 
and one may also expect topological superconductor-like 
states may be described by similar models in the main 
text. However, the exactly model is unable to construct 
because the superconducting pairing is a many-body 
effect and we have to rely on some mean- field treatment. 
A straight way to obtain a topological superconductor- 
like state is to introduced attractive local interaction 
Hu = -C^SjC-t'^itcI^c^ into Eq.lgni). According to 
Ref.[53[, assuming a s-wave paring and solving the 
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corresponding mean-field self-consistent equation for 
the superconducting order parameter, one can obtain 
an effective p + ip chiral superconducting phase with 
non-zero Chern number Ci. More details can be found 
in Ref . [50| . Then, combing the electron part with the Z2 
gauge-field, we expect there exists a transition from the 
chiral superconducting phase to the Z2 fractionalized 
chiral superconductor. 



charge/spin Hall conductance) are identical with single 
particle excitation different. An extra element for the 
present orthogonal Chern/topological insulators, which 
is not noticed in original studies^^^— comes from their 
edge-states for physical electrons as shown in the cases 
for the square and honeycomb lattices. 



D. Relation to the fractional Chern insulators 

It is also interesting to discuss the relation of Z2 frac- 
tionalized Chern insulators in our paper to the fractional 
Chern insulators in Refs.f32l-t36|. We note that in their 
model, the filling of electrons is fractional, thus the sys- 
tem is metallic without introducing interactions while our 
case is insulating since half-filling is explicitly assumed. 
To achieve the fractional quantum Hall-like states, the 
noninteracting part of the Hamiltonian is finely tuned 
into nearly-dispersionless (flat) as the case for the flat 
Landau levels in the classic fractional quantum Hall ef- 
fect. Then, including the effect of certain proper interac- 
tions leads to fractional quantum Hall-like states found in 
their sophisticated numerical calculation. For our mod- 
els, the Landau levels-like element is not involved, thus 
we do not require the flat band condition. The fraction- 
alized insulating state in our models will have Z2 topo- 
logical order in contrast to the m-iold (to — l/v with 
ly being the particular fractional filling) topological de- 
generacy in the fractional quantum Hall-like states. In 
other words, on a torus, Z2 fractionalized Chern insu- 
lators discussed in the main text have four degenerate 
ground-states while m degenerate ground-states appear 
in the fractional Chern insulators. Therefore, these two 
kinds of fractionalized states are quite different in their 
nature and it seems still a challenge to construct even 
artificial exactly soluble models for such fractionalized 
Chern insulators. 



E. Relation to the orthogonal metals 

The orthogonal metals are metallic in their original 
definition but it can also extend to gapped phases 4^"— 
When such extension is made, one expects a fractional- 
ized state and a non-fractionalized counterpart. These 
two phases should have the same bulk thermal and 
transport properties. From the discussion in the pre- 
vious sections, we realize that the Z2 fractionalized 
Chern/topological insulators in exactly soluble correlated 
models can correspond to the fractionalized state and 
the usual Chern/topological insulators relate to the non- 
fractionalized ones in orthogonal metal-like states. In 
the sense of orthogonal metals, the Z2 fractionalized 
Chern/topological insulators can also be dubbed as or- 
thogonal Chern/topological insulators since bulk thermal 
(gapped excitations) and transport properties (quantized 



VI. CONCLUSION 

In summary, we have propose certain exactly 
soluble models which support Z2 fractionalized 
Chern/topological insulators besides the usual 
Chern/topological insulating states. The bulk be- 
haviors of physical electrons are similar in these two 
states but the edge-states of physical electrons have 
rather different behaviors, which provides a definite 
signature to identify the fractionalized states from the 
non-fractionalized ones. Besides, the transition from the 
usual Chern insulator to the Z2 fractionalized Chern 
insulator is found to fall into the usual 3D Ising universal 
class. 

Moreover, we have inspected relations to the QSH* in 
RefUTI and have made an extension to the case of the 
square lattice. For the specific case of the square lat- 
tice, we expect that a chiral topological superconducting 
phase and its Z2 fractionalized version may appear when 
the attractive local interaction is introduced. We also 
provide intimate link of our exactly soluble models to 
various more real lattice models without intrinsic degree 
of freedom for gauge fields. In terms of the Z2 slave- 
spin mean-field approximation ; ^ ^ 1 ^^ we have suggested 
that similar fractionalized states found in our present 
paper may appear in those lattice models. The rela- 
tion of Z2 fractionalized Chern insulators found in the 
main text to the fractional Chern insulator Refs. |33 - [3a | 
is inspected, which implies that these two kinds of frac- 
tionalized states are quite distinct in their nature. We 
also provide some discussion on the relations to orthog- 
onal metalsi^"— and fractional topological insulatorsJ^ 
The present work may be helpful for further studies on 
the fractional Chern/topological insulator and the re- 
lated novel strongly correlated quantum phases. 
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Appendix A: The standard Zi lattice gauge theory 

model on the honeycomb lattice and its topological 

order 



The standard Z-i lattice gauge theory model on the 
honeycomb lattice is defined by the following Hamilto- 



B. 



-^-n'jY. n ^n 



(Al) 



(u> 



i j^hexagon 



where the first term describes the fiuctuation induced by 
afj while the second one denotes the potential gained by 
six Ising field af in the same hexagon. One can define 

the Z2 gauge transformation operator Gi — ni=i+(5„ ^^j 

and it is easy to check that [Gi,Hz2] = 0. Thus, the 
model Hz2 has the wanted Z2 gauge invariance and any 
physical states should be invariant under the operation 
of Gi. There are two limit situations which are easy to 
analyze. First, if J ^ /i, then the ground-state should 
have Ylj^hexagon ^ij = ^ ^'^^' ^^^ hcxagous and this is just 
the deconfined state for the Z2 gauge theory. However, 
when J <^ h, one may set all af = 1 to reach the ground- 
state and the obtained state is the so-called confined state 
for the Z2 gauge theory. 

Then, we briefly discuss the topological order of the 
Z2 lattice gauge theory in its deconfined state.^ Consid- 
ering a cylindrical geometry and assuming the we work 
in the deconfined phase, the ground-state is easily ob- 
tained by setting all of afj = 1 and this corresponds to 
the case without threading a Z2 vortex (vison) (A Z2 vor- 
tex is created by setting trf- = — 1 in certain link while 
keeping others unchanged.) into the hole of the cylinder. 
If we thread a Z2 vortex into the hole of the cylinder, the 
resulting state is also the ground-state. Thus, the decon- 
fined phase of the Z2 gauge theory has two degenerate 
ground states. Using the same argument, the the decon- 
fined phase of the Z2 gauge theory on a torus will have 
four degenerate ground states corresponding to the Z2 
vortex threading or not threading each of the two holes. 

The benefit of the above model is that there exists 
an exact dual transformation which transforms the Z2 
lattice gauge theory model into a quantum Ising model 
on the honeycomb lattice. The dual transformation is 
defined as 



n 



j^hexagon 



<J,,j = T, Tj . 



(A2) 



(A3) 



Then, the original Hamiltonian is transformed into the 
following one 



(ij) i 



(A4) 



One may note that this is just the quantum transverse 
Ising model defined on the honeycomb lattice. Reader 



may care whether some constraints are needed for the 
transformed Hamiltonian H. It is crucial to see that the 
Z2 gauge transformation operator Gi — Y[j=i+Sa ^ij ~ 
Y[-i=i+s ^j ^'^'^ ^^^^ should be taken as a constraint in the 
dual quantum Ising model. Therefore, in general, the re- 
sulting quantum transverse Ising model is not free. If we 
simply neglect the effect of the constraint, the quantum 
transverse Ising model will have two phase characterized 
by the Landau local order parameter (ff ) = (the mag- 
netically disordered phase) and (ff ) ^ (the ferromag- 
netic phase). There also exists a second-order quantum 
phase transition between those two phases when t/ J is 
finely tuned into the critical value. The critical expo- 
nents of this quantum model is identical to the classical 
Ising model in three space dimension, [a = 0.11, /3 = 
0.32, 7 = 1.24, 5 = 4.9, v = 0.63, 77 = 0.04, z = 1.) 



Appendix B: Path integral and effective theory for 
the quantum transverse Ising model 

The quantum Ising model in transverse field is defied 
9 



2 
as^ 



H, 






-^E- 



(Bl) 



where a ferromagnetic coupling /i > is assumed and J 
represents the the transverse external field. 

At first glance, one may directly use the coherent state 
of spin operators in constructing the path integral rep- 
resentation, (One can find a brief but useful introduc- 
tion to this issue in Ref. [2]) however, this will lead to 
an extra topological Berry phase term and is not easy 
to utilize practically. An alterative approach is to use 
the eigenstates of spin operator r^ or r^ as the basis for 
calculation.— One will see this approach is free of the 
topological Berry phase term and give rise to a rather 
simple formalism. Therefore, to construct a useful path 
integral representation, we will follow Ref. [60| . 

First of all, we consider the orthor-normal basis of Ns- 
Ising spins as 



\a) = \(Ji)\a2)\(T2) ■ ■ ■ \(Jn) 
with ai = ±1 and define 



\o) 



Ti |0-) = \(n)\(T2)\CT3.) 



o-,: o- 



(Ji) ■ ■ ■ \ctn) 



(B2) 



(B3) 



(B4) 



Then the partition function Z = Tr{e ^^) can be rep- 
resented as 



N 



{(T} = ±ln=l 

where eN=/3. The calculation of {c7{7i + l)\e'-^^^''i \a{n)) 
is straightforward by exponentiating the rf matrix and 
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one gets 

(a(n + l)|e''^^.^n^W)= ^(e"^ + e-^'V,(n)a,(n + 1)), 

where a = i[lncosh(eJ) — Insinh(eJ)] and h = 
■i[lncosh(eJ) + lnsinh(eJ)]. Therefore, the resuhing path 
integral formahsm for the quantum Ising model in trans- 
verse field is 

N 
{cr}=±ln=l 

(B6) 

Further, if one assumes the model is defined in a hyper- 
cubic lattice in space dimension of d, an effective theory 
can be derived as 

Z= /i?0<5(02_l)e-/''-'^''-*[(^^^)'+-'(^«'l, (B7) 



where 5^- = (-t-) 2 with oq being the lattice constant 



and c 



2s ^a, 

I, d 

2 „ 'iOq 



Moreover, in the effective theory, </> cor- 



responds to T^ while T^ gives the kinetic energy term 
in imaginary time. Then, the standard 0^ theory is ob- 
tained by relaxing the hard constraint (fp' = \ while in- 
troducing a potential energy term. 



D(f)e 



' f dTd'*£c[(a^0)^+c2(V0)^+r02+„04j 



(B8) 



where r, u are effective parameters depending on micro- 
scopic details. 

For the case of the honeycomb lattice, since only the 
modes near the minimum energy region of the band are 



involved in low-energy limit, the quantum Ising model on 
the honeycomb lattice can also be described in terms of 
the above ip'^ theory in spite of its bipartite feature. 



Appendix C: Derivation of Chern-Simon action 

Here, we would like to give a brief derivation of the ef- 
fective Chern-Simon action Eq. (jH]) from the Dirac action 
Eq. ([5]). First, the Dirac action is written as 

S= d'^xd.T^liPaaht.id^, - ieA^) + m)'il;aa]-{Cl) 

Then, integrating out Dirac fermions one obtains 

Seff = TV In ZJet [7^(9^ - ieA^,) + m] 

= 7VTrln[7^(a^-ieA^)+m] 

= 7VTr[ln[7^a^ + m] + ln[l - ze(7^a^ + m)-^-f^,A^]] 

c^ N J ^A,U,Aq)A. (C2) 

wtierelV(<7j - ^- j -(2^^r[ ^i_^_^i 7^ ^'i+(^k+qyi ^mJ 

^ / (0 L4fc^ „.+(!.+,). ][-»"-tgArr(7^7.7A)] + ... ~ 
em^/ii/A^^ is calculated at one-loop level. We 

have also used the identity Tr{'^^'y^"f\) — 2ie'^''^ 
with 70 = Tz, 71 = Tx and 72 = Ty while 



\/q^+4m^ 



for 



r d^k r 1 1 1 _ ^^^ 

J JIF)^ \-rnF+k^ m^+ik+qp' i 4Tr\q\ — 87r|m| 

l^l <C \m\. Therefore, the effective Chern-Simon action 
Eq. ([5]) is obtained as 



Seff- 



- — T^NAn — ; — r' 
(27r)3 ^87r|m| 



/11/A 



q\A^ 



= / d^xdrlNe^ 



87r|?n 



■e'^-'^A^d.A^] 



(C3) 
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